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SUMMARY 


An analytical study of hydrogen-air kinetics has been performed as part 
of the scramjet engine research efforts at Langley Research Center. Calcula- 
tions were made over a range of pressure from 0.2 to 4.0 atm, temperatures 
from 850 to 2000 K, and mixture equivalence ratios from 0.2 to 2.0. The 
finite-rate chemistry model included 60 reactions in 20 species of the 

calculations also included an assessment of the effect 
of small amounts of the chemicals ^ 20 , ^2^2' ^3 t-he, initial mix- 

ture on ignition and reaction times as well as an assessment of the effect of 
the variation of the third-body efficiency of H 2 O relative to N 2 in certain 
key reactions on reaction time. 

The results indicate that for mixture equivalence ratios between 0.5 and 
1.7, ignition times are nearly constant; however, the presence of H 2 O and NO 
can have significant effects on ignition times, depending on the mixture tem- 
perature. Reaction time is dominantly influenced by pressure but is nearly 
independent of initial temperature, equivalence ratio, and the addition of 
chemicals. Effects of kinetics on reaction at supersonic combustor conditions 
are discussed, and recommendations are made for further experimental studies 
of rate constants for key reactions. 


INTRODUCTION 

The study of the fluid dynamics and combustion processes of hydrogen- 
fueled supersonic combustion ramjets (scramjets) has been a major part of the 
hypersonic propulsion research program at Langley Research Center. Most of 
these studies have been, of necessity, experimental investigations of the 
overall effects of fuel injection, mixing, and reaction on the flow in sub- 
scale engine components. This dependence on empirical results is due to the 
complexity of the flow around fuel injectors in three-dimensional geometries, 
which are not readily treated analytically. Most numerical solutions are 
restricted to two- or three-dimensional parabolic flow with equilibrium or 
simple finite-rate chemistry models of the H 2 '“air system. Such numerical 
solution schemes are applicable only in the parabolic flow region well down- 
stream of the disturbance caused by the transverse fuel injection. Scramjet 
combustors rely on transverse injection of the fuel in order to achieve rapid 
mixing and reaction. To begin the calculation therefore requires a prior 
knowledge of the extent of fuel mixing, ignition, and reaction. 

All this information, particularly the ignition-reaction state, may not 
even be available from experimental data because of the difficulty of freezing 
the reaction in a sample extracted by a probe from the high-velocity, high- 
temperature combustor flow. In spite of these problems, however, a sufficient 
data base has been established to define a scramjet engine concept and to per- 
mit fabrication of subscale engine models with integrated inlet, combustor, 
and nozzle components. The current scramjet is designed for operation over a 



flight Mach number range from 4 to ?• The corresponding stagnation tempera- 
tures of the captured air are 900 K and 2200 K, respectively, which should be 
sufficient for autoignition of hydrogen in the separated flow regions near 
fuel injectors. In recent ground tests of subscale engine models (ref. 1), 
however, problems were encountered in obtaining ignition and sustaining reac- 
tion at test conditions for which ignition and reaction might have been 
expected. The need to better understand the chemical mechanism of the igni- 
tion and reaction of H 2 ”air mixtures at conditions typical of a scramjet com- 
bustor led to this study. 

The present analytical study uses a computer program (refs. 2 and 3) for 
solving chemical kinetics problems in one-dimensional flow, with an H 2 -O 2 -N 2 
finite-rate chemistry model consisting of 60 reactions involving 20 species. 
The calculations covered a range of pressures from 0.2 to 4 atm, a range of 
temperatures from 850 K to 2000 K, and H 2 “air mixtures between 0.2 and 2.0 
times the stoichiometric value. In addition, the effects of contaminants such 
as NO^ and H 2 O, which may be present in the high-temperature test gases used 
in ground test facilities, were investigated. Significant concentrations of 
NO^ molecules are produced in facilities which use an electric arc to heat the 
test air; H 2 O is produced in combustion-heated facilities as a result of burn- 
ing hydrogen in air enriched with oxygen so that the airlike products have the 
proper oxygen mole fraction. (In the latter facilities, the test gas product 
is sometimes referred to as "vitiated" air. ) The effects on ignition times of 
potential ignition aids such as H 2 O 2 and O^, which may be added to ground test 
facilities, were also considered. In addition, is naturally encountered in 
high-altitude (>30 km) flight. Because the computer program is one- 
dimensional, the effects of fuel injection and mixing cannot be included. 
However, the present finite-rate calculations of premixed H 2 ~air flows are 
considered important in predicting ignition and reaction in the recirculating 
zones near the fuel injectors. In analysis of scramjet combustors, where the 
reaction may be limited by the mixing, these results for premixed reactants 
provide a lower limit of the times required for reaction. 

Many other analytical studies of the chemical kinetics of the H 2 "air 
system, including vitiated air, have been conducted (e.g., refs. 4 to 19). 
However, many of the problems encountered in scramjet testing have not been 
addressed explicitly. These include ignition and reaction at low tempera- 
tures; varying amounts of water in the vitiated-air test gas; and realistic 
equilibrium concentrations of oxygen and hydrogen atoms, hydroxyl free radi- 
cals, and oxides of nitrogen in the test stream. For example, in the vitiated 
test gas of a ground test facility at conditions representative of Mach 4 and 
Mach 7 flight, water mass fractions are about 0.05 and 0.23, respectively. In 
an arc-heated facility, NO^ mass fractions may range between 0.01 and 0.03. 

In some of these previous calculations, the chemistry models were too simple 
to be applied at low-temperature or high-pressure conditions. In some refer- 
ences, the reaction-rate constants were either estimated, because of a lack of 
data, or they were approximated from early data that precede the survey by 
Baulch et al. in references 20 and 21. This survey (refs. 20 and 21) along 
with more recent data (ref. 15) was used to update the rate constants of the 
H 2 “air mechanism. 
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The purpose of the present one-dimensional calculations is to define the 
conditions at which chemical kinetics are important in future ground tests of 
subscale scramjet-engine models and in analyses of the combustor flow using 
more sophisticated computer codes. The present calculations include (1) the 
effects of initial temperature, initial pressure, and the equivalence ratio 
(degree of stoichiometry) of the H 2 "air mixture on the ignition and reaction; 
(2) the effects of contaminants such as H 2 O, OH, H, and O present in an equi- 
librium vitiated test gas; ( 3 ) the effects of NO^ on ignition at low tempera- 
ture; ( 4 ) the effects of H 2 O 2 and O 3 used to promote ignition in ground test 
facilities and the effect of atmospheric O 3 that may be encountered in scram- 
jet flight; and ( 5 ) the sensitivity of ignition and reaction to the reaction 

rates of dominant reactions. 
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SYMBOLS 


area, m 



(19) ), 


species production function (eq. ( 11 )) 


activation energy in rate-constant equation (eq. { 19) ) 


enthalpy production function (eq. ( 12 )) 

molar heat capacity of species i, J-mol 

specific heat at constant pressure, J-g -K 

local and stoichiometric hydrogen-to-oxygen mass ratios 

molar enthalpy of species i, J-mol 

total-mixture static enthalpy per unit mass, J-g 

total-mixture energy per unit mass, J-g ^ 

equilibrium constant 

backward reaction-rate constant 

forward reaction-rate constant 

number of chemical reactions 


Mach nvimber; third-body species in certain reactions 
third-body efficiency factor defined by equation (25) 
molecular weight of species i, g-mol ^ 
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mass flow rate, g-s*' 

third-body efficiency factor of species i in reaction i 
(see eq. (26) ) 

number of species in gas mixture 

temperature exponent in rate-constant equation (eq. (19)) 
pressure, atm (1 atm = 101.325 kPa) 

universal gas constant, 1.987 cal-mol”^-K"^ (1 cal = 4.184 J) 
species i (eq. (18)) 
temperature, K 
time , s 

velocity, m-s"^ 

net species production rate (eq. (21)) 
mixture molecular weight, g— mol""'^ 
mass fraction of species i 
ratio of reaction times 
specif ic-heat ratio 
ratio of ignition times 

forward stoichiometric coefficient in general reaction 
(eq. (18)) 

reverse stoichiometric coefficient in general reaction 
(eq. (18)) 

density, kg-m~^ 

concentration of species i, moles per unit mass of mixture, 
mol-g ' 

ignition time, s 
reaction time, s 
total reaction time, s 
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equivalence ratio, f/f* 

net conversion rate in kinetics equations 

net rate of formation of ith species in reaction j 





Subscripts : 

eq 

i 

j 

o 

r 

Abbreviations : 


equilibrium 
species (i = 1/ 2, • 

reaction (j = 1, 2, 

initial conditions 
reference conditions 


. N) 

i) 


EFS 

ppmw 


equilibrium free stream 
parts per million by weight 


ANALYSIS 

Differential Conservation Equations 

The computer program used for this analysis of a complex reacting gas 
mixture flowing through an arbitrarily assigned area employs the conservation 
equations for one-dimensional steady-state flow. The flow is assumed to be 
isobaric, premixed, inviscid, and adiabatic. Additional information about 
these equations and the numerical integration procedure used in the computer 
code may be found in references 2 and 3. The equations employed are 

Global mass conservation 


pAu = m 


( 1 ) 


Continuity for each of the N species in the mixture 



W. Au 
1 


(i = 1, 2, . . .r N) (2) 


where O is the concentration of species i, in moles per unit mass of the 
mixture, ""and W. is the species production rate. Expanding^ the differentia- 
tion of equation (2), using equation (1), and noting that dm/dt - 0 reduces 
the equation to 


1 / 2 , . 
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Conservation of momentum 
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u dt 


= 0 
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Conservation of energy 
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(5) 


where xs the total constant energy per unit mass of the gas and the mix 

ture molecular weight is defined as 


W 

m 




( 6 ) 


Equation of state 


PRT 
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(7) 


By use of equations (1) to (7), the following differential equations can be 
derived, with time t the independent variable: 
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Other parameters used in these equations are defined as follows: 
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(17) 


Equations (3) and (8) to (10) form a system of N + 3 equations in N + 3 
unknowns. This equation system is solved by numerically integrating with time 
for a specified area distribution. The equation set for a specified pressure 
distribution and for the passage of a normal shock through the gas is given in 
reference 2. 


Chemical Reaction Equations 

Each chemical reaction can be written in the general form 


v-r ^ ■ 

N D N 

y) V. . s. + yi V,' . s. 

i=i 1 1 


(18) 


where represents the ith species. The forward reaction rate is given by 


kj = AjT ^ exp(-Bj/RT) 


(19) 


The backward reaction rate is determined from k- and the equilibrium con 
stant Kj ; ^ 


k . 

” K. 
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( 20 ) 


The rate of formation of each species which appears in equation (3) is 

defined as 
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ij 


( 21 ) 
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where 


a)=p^(v'-v)\(;. (22) 

ij ij ID D 


is the net rate of formation of species i by reaction j. The parameter 

is the net reaction conversion rate introduced in reference 2 as an indi- 
cation of the amount of change due to each reaction that occurs in the mix- 
ture. The form of depends on the type of reaction. Three types of 

reversible chemical reactions are considered in the computer program: 


Type 1 - Bimolecular shuffle reaction 


k . 

^ ^ ^4 


with 


D 




(23) 


Type 2 - Three-body recombination 


k. 

M + S^ + S^ ^ 3 "*” ^ 

(where M can be any species present) with 


= k .M . 
D D D 


pa a 
. 1 2 



(24) 


In this equation, Mj is the third-body efficiency factor for the jth reac- 
tion defined as 
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m. . 0. 
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( 25 ) 


M. 
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N 

E 
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and is the third-body efficiency factor for species i in reaction 

j. This factor is a correction to the preexponential factor A. in the rate 
equation (19) and can be written as ^ 


in. 

^f3 


(A.) 

'3 s, 


(A. ) 

3 r 


(26) 


where (Aj)j~ is the value of Aj for a reference species as third body. For 
the calculations presented in this paper, all third-body efficiencies are 
referenced to N^. Except as listed at the bottom of table I, all third-body 
efficiencies are 1.0. 

Type 3 - Two-body dissociation 


M + 



+ s^ + M 
2 


with 


if, = k .M. 
3 3 3 



-r) 


(27) 


For some problems the heat release is an important consideration. In these 
situations a useful quantity is the net energy conversion rate for the jth 
reaction defined as 


(ifj . = ip.(AH 
^3 3 


298 


j 


(28) 


where (AH 29 q)^ is the molar heat of reaction at 298 K for the jth reaction 
proceeding in the forward direction, from left to right. 
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Computer Program 


The computer program for solving the one— dimensional flow with finite- 
rate chemical reactions was first presented in reference 2. The system of 
conservation equations and chemical reactions form a nonl inear ^ coupled set of 
differential equations, which are solved by the Adams implicit method 
(ref* 22). A later modification of the computer code (ref. 3) added the 
stiffly stable, linear, multistep method of Gear (ref. 23). The incorporation 
of the Gear method into the program provides a solution technique that is a 
highly efficient means of altering step size to obtain a stable rapid solution 
to the large system of stiff equations. This is of particular advantage when 
the number of chemical species is greater than 15. 

Input to the computer program is user oriented. This allows for easy 
modification to the chemical reaction system, reaction-rate parameters, ther- 
modynamic properties (ref. 24) and general problem specification. Additional 
details of the input specifications are given in reference 3. 


Scope of Present Calculations 

The objective of the present study is to parametrically examine the igni- 
tion and reaction times of H2~3,ir mixtures at conditions representative of 
scram jet combustor flows. The ignition time, also known as ignition delay or 
induction time T. , is the time required for the rapid buildup of the free 
radicals and atoms^with very little heat release. Three definitions of igni- 
tion time have been used in previous research: (1) the time for the concen- 

tration of OH to reach 5 x 10“^ mol/cm^, (2) the time at the end of the con- 
stant exponential growth of OH mass fraction, and (3) the time taken to reach 
a temperature at ignition T^^ of 

T. = 0.05(T^ - T^) + (29) 

ig eq o o 


The ignition time given by definition (1) was used only for comparison with 
the experimental data of reference 15. It was observed that the ignition 
times given by definitions (2) and (3) are approximately the same over a wide 
range of pressure and temperature. The ignition— time results from these 
calculations are therefore presented using equation (29). The total reaction 

time is defined as the time taken to reach a temperature T,pj^ where 

TR 


T 

^TR 


0.95(T 


eq 




(30) 


Reaction time is then defined as 


T 
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T 

TR 


ig 


1 1 


(31) 



Chemistry model .- The chemical reaction model, given in table I, consists 
of 60 (£ = 60) reactions involving 20 (N = 20) species: 


02 

N2 

«2 

0 

»2° 

OH 

HO2 

NO 

CO 

CO2 

HNO 

HNO, 

H2O2 

O3 

HCO 

Ar 


N H 

NO2 

HNO3 


Reaction rate data were taken from references 14, 20, 21, and 25 to 32. These 
data are some of the latest available for the hydrogen oxidation mechanism and 
provide improvements to the reaction-rate parameters and third-body efficiency 
factors for certain key reactions such as 


H + O + M HO + M 
NO + HO^ -»• OH + NO^ 

and 

H + OH + M(H 20 ) > H^O + M(H 20 ) 


The impact of the reaction rates and third-body efficiencies of these reac- 
tions on ignition and reaction is discussed in the section "Reaction Time." 

Parameter variation.- The principal parameters of these calculations are 
the initial static pressure and temperature and the hydrogen-air mixture com- 
position. Static pressures were varied between 0.2 and 4.0 atm; static tem- 
peratures between 850 and 2000 K. The amount of hydrogen in the initial 
mixture was varied corresponding to equivalence ratios between 0.2 and 2.0. 
The equivalence ratio <j> is an indication of the relative stoichiometry of 
the mixture based on the global reaction for the oxidation of hydrogen: 


2 H 2 + 2n^0 


For this reaction the stoichiometric mass ratio of hydrogen to oxygen is 


(moles H^) (MW) 

^ _ 2(2.016) _ 

(moles 02 )(MW) 1(32) ~ 0.126 ( 32 ) 


For an arbitrary mixture of hydrogen and air (or other species) with a mass 
ratio of hydrogen to oxygen f, the equivalence ratio is defined as 
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( 33 ) 



Most of the calculations were made for <J) = 1.0. 

Chemical additives and contaminants In addition to the variation in the 
state properties, the composition of the air was varied to study the effects 
of contaminants {H 2 O, present in ground test facilities and chemical 

additives (H 2 O 2 , O^) used as potential ignition aids in ground tests of sub- 
scale scramjet engines. In order to simulate the combustor conditions of a 
scram jet at Mach 4 to Mach 7 flight, it is necessary to heat the air in ground 
test facilities. One way to produce a high-temperature test gas is to burn 
hydrogen in oxygen-enriched air. The test gas produced by such a combustion 
heater facility (ref. 33) is often referred to as vitiated air, since it has 
the same oxygen mole fraction as air but is contaminated by a significant con- 
centration of water vapor. The amount of H 2 O present in the vitiated test gas 
is directly related to the desired stagnation temperature (or flight Mach num- 
ber) that is simulated. For this analyis, the following values of water mass 
fractions are assumed to form a consistent set with each of the initial 
temperatures : 

T ,K 910 1000 1100 1200 >1300 

o' 

OL, ^ 0.097 0.125 0.145 0.20 0.245 


Whereas previous analyses have accounted for vitiation effects by arbitrarily 
varying the concentraton of H 2 O and OH free radicals at the same initial tem- 
perature, the consistent pairs of water mass fractions and initial tempera- 
tures used in the calculations here provide a more realistic means of assess- 
ing the effects of vitiated air on ignition and reaction. 

Another scramjet ground test facility (ref. 1) uses an electric arc to 
heat the air and produces significant amounts of nitrogen oxides NO and N02* 

It has been estimated that this facility produces total NO^ mass fractions of 
about 0.01 and 0.03 at conditions simulating Mach 4 and Mach 7 flight, respec- 
tively. Calculations were made to ascertain the effects of NO and NO 2 mass 
fractions of 0.01 on the ignition times of the H 2 *“air mixtures. 

In addition, the effects of free radicals that may be present in the test 
gas as a result of dissociation at the high temperatures were examined. The 
amount of each free radical present was determined by assuming that the air or 
vitiated air of the mainstream test gas was in chemical equilibrium at the 
initial pressure and temperature. Thus, for calculations of H 2 in equilibrium 
air, small amounts of O, NO, and NO 2 were included in the initial mixture. 

For H 2 in equilibrium vitiated air, small amounts of O, H, OH, HO 2 / NO, and 
NO 2 were included in the initial mixture. In the following discussion, the 
results obtained with equilibrium free-stre 2 um (EFS) mixtures are referred to 
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as H 2 -air (EPS) and H 2 ~vitiated air (EPS) to denote hydrogen mixed with equi- 
libria air and equilibrium vitiated air, respectively. In the vitiated-air 
facility, the test gas would likely be in equilibrium at the local temperature 
and pressure at the combustor entrance. However, the arc-heater facility 
would produce high levels of (mass fractions of 0.01 to 0.03) that would 

freeze as the flow expanded through the facility nozzle, so that a true equi- 
librium would not be obtained. 

The effect of small amounts (mass fractions of 0.01) of ^ 2^2 
ignition was examined to assess their potential as ignition aids in gro^d 
tests. Calculations were also made with an O 3 concentration of 10 pp«nw to 
approximate the atmospheric levels encountered in flight at an altitude of 
30 km. Maximum 0^ concentration in the atmosphere is about 13.5 ppmw at 36 km 
(ref. 34). 


RESULTS AND DISCUSSION 

The one--dimensional kinetics computer code (refs. 2 and 3) previously 
described was used to calculate ignition and reaction times for the range of 
parameters indicated. Typical output of the computer code is presented in 
figure 1 in the form of the time variation of species mass fractions through 
the ignition and reaction zones. Results are given for a stoichiometric 
air mixture at an initial temperature of 1000 K and at initial pressures of 
0.5 and 2.0 atm in figures 1(a) and 1(b), respectively, and for a stoichiomet- 
ric mixture of H 2 in equilibrium vitiated air at T^ = 1000 K and p = 1 atm 
in figure 1(c). The times of ignition and total reaction are noted on each 
figure by and respectively. The effect of increased pressure to 

shorten reaction time can be seen by the much smaller difference between x 

^ig figure 1 (b), at 2.0 atm, when compared with the difference indi- 

cated in figure 1 (a), at 0.5 atm. 


Ignition Times 


The effects of initial pressure and temperature on ignition time of a 
stoichiometric H^-air mixture are presented in figures 2 (a) and 2 (b), respec- 
tively. Because of the widespread use of ignition times from reference 4, 

^ig ^^sults as given by reference 4 are plotted for comparison as dashed 
lines in figure 2(a). In contrast to the constant slope with pressure of the 
reference 4 ignition times, the present results indicate a drastic nonlinear- 
ity in the pressure dependence, particularly at values of T below 1100 K. 

At T^ above about 1200 K, the present results indicate a slight retardation 
of ignition times relative to those given in reference 4. The constant- 
pressure curves of x^^ as a function of the reciprocal temperature in fig- 
ure 2(b) indicate very long ignition times as T decreases. These effects 
produced by the present chemical mechanism are primarily due to the formation 
of the HO 2 molecule through the reaction H + O 2 + N 2 > HO 2 + N 2 , which was not 
included in the kinetics model of reference 4 . This reaction is a chain 
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terminating reaction which depletes H atoms* Because it is a three-body reac- 
tion, it should become more important at higher pressures. Since the reaction 
rate is weakly dependent on temperature, this reaction is important at low 
temperatures, where the chain carrier reactions are slower. 

The sensitivity of ignition times to the fuel equivalence ratio was 
checked for and H 2 “Vitiated air at 1200 K and 1.0 atm. Results of 

these calculations are given in figure 3 along with previous results for H 2 ~ 
air from references 8, 16, and 19. Also included in figure 3 are data for 
1111 K and 1389 K taken from reference 18 in the format presented in refer- 
ence 17. The present ignition times agree reasonably well with the referenced 
results, a spread of about 20 percent in occurring for (f) between 0.5 

and 1.7. In comparison with the data, the present predictions of ignition 
time at 1200 K show the same trend with increasing and fall between the 

data for 1111 K and 1389 K. For both H 2 “air and H 2 “Vitiated air, ignition 
times predicted by the present chemistry model exhibit only a ±10 percent 
variation for (f) in the range from 0.5 to 1.7, although the ignition times 
for H 2 “vitiated air are longer than for H 2 “air by about 20 percent. This 
result supports the conclusion of reference 4 that ^±g constant for 

between 0.4 and 2.0. 

In order to test the correctness of the ignition times predicted by the 
present chemistry model, the theoretical results have been compared with stoi- 
chiometric H 2 “air shock-tube data from reference 15. These comparisons are 
presented in figure 4(a) for initial pressures of 0.5 and 2.0 atm and in fig- 
ure 4(b) for initial pressures of 0.27 and 1.0 atm. The agreement between 
theory and data is quite good, except for p = 0.27 atm in figure 4(b). The 
authors of reference 15 reported that these low-pressure data were the least 
reliable because of boundary-layer growth in the shock tube. 

A cross plot of these reference 15 data at = 1000 K is given in fig- 

ure 5 with theoretical predictions from references 4 and 14 and the present 
calculations for comparison. The present results agree well with the data 
except at the low-pressure (p = 0.27 atm) datum. The theory of reference 14 
includes a sophisticated chemistry model with 25 reactions. It includes 
hydroperoxyl (HO 2 ) and hydrogen peroxide (H 2 O 2 ) molecules and shows the proper 
trend, although it does not agree with the data. Since the chemistry model of 

reference 4 did not include the HO 2 reactions, it would not be expected to 

agree with the data at high pressures (>1.0 atm), where the HO 2 molecule 

becomes important. The linear relation between and p of reference 4 

does, however, agree with the data in the narrow range of pressures between 
0.5 and 1.0 atm at T^ = 1000 K. 

Effect of vitiated air .- Some scram jet ground test facilities (ref. 33) 
rely on the combustion of H 2 to produce the desired high-temperature test 
gas. Because the resulting vitiated-air test gas is contaminated with high 
concentrations of H 2 O, the effect of this contaminant and its dissociation 
products on ignition times has been examined. A comparison of ignition times 
for stoichiometric H 2 -air (EFS) and H 2 “vitiated air (EFS) is presented in 
figure 6. In figure 6(a), in which the effect of T^ on the variation of 
X. with p is presented, it can be seen that at T^ less than 1100 K and 
p^^greater than 0.2 atm, the presence of H 2 O in the vitiated test gas retards 
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ignition, particularly at the lower value of T^. At temperatures above 
1200 K, the effect of the vitiation is to enhance ignition. The reason for 
the longer ignition time at < 1100 K is that in a vitiated-air stream, 
the production of HO 2 is accelerated by the reaction H + O 2 + H 2 O HO 2 + H 2 O 
because the third-body efficiency of H 2 O may be as much as 13 times that of 
N 2 . Although some OH radicals are present in the free stream, the foregoing 
reaction tends to deplete H atoms , an important chain carrier in the ignition 
region, and cause an increase in ignition times. 


At temperatures above 1200 K, ignition is enhanced by the dissociation of 
H 2 O in the vitiated free streeim and the resulting production of OH radicals, 
whose concentration (10“^ to 10"^ by mass) is significantly greater than the 


mass concentration of O and H atoms. The curves of x. versus reciprocal 

in figure 6(b) also illustrate the 


temperature (1000/T^) at constant p 


effect of a vitiated free stream on ignition times. Other calculations 
(refs. 7 and 16) exhibit similar effects at low T^ over a limited pressure 
range. Therefore, at low initial temperatures, the ignition characteristics 
observed in ground tests of scramjet combustors using a vitiated stream are 
longer than those for tests with air at similar conditions. At the higher 
initial temperatures, the ignition times are shorter with a vitiated stream. 


Effect of equilibrium free stream .- Figure 6 compares the difference 
between air and vitiated-air free streams when both are in chemical equi- 
librium. The effect on x^^ of the free radicals produced by the dissocia- 
tion of H 2 O in the vitiated air at chemical equilibrium can be seen from the 
curves in figiore 7. This comparison is made for <(> = 1.0 at T = 1300 K. 

For H 2 “air (solid lines), the presence of O and NO in an EFS (otQ = 1.3 x 10 "®, 
= 4 X 10 ) has little effect on the ignition times. For the H 2 ~vitiated 

air (dashed lines), in which ^ = 0.245, the assumption of EFS at 

” 2 ° 

Tq = 1300 K ^^roduces OH, 0, H, and NO («qjj = 1*5 x 10~^, = 6 x 10~®, 

otjj = 1 X 10 , = 3.5 X 10 which reduce ignition times by about 

40 percent. Note that for the vitiated air in which the H 2 O is not dissoci- 
ated (not in equilibrium) , the values of ^iq these conditions are larger 

than for air. These comparisons suggest that care must be exercised when 
applying theoretical results to insure that the effects of vitiation are 
properly accounted for. The investigation reported in reference 35 emphasized 
the necessity of considering free radicals when comparing theory with combus- 
tion data for parallel H 2 injection into a vitiated test stream. Reference 8 

also indicated a sensitivity of T . to the amounts of H, O, and OH. 

ig ' 

Effect of HpO. - The presence of H 2 O in vitiated air has been shown to 
strongly affect calculated ignition times, particularly at initial tempera- 
tures of 1000 K or less. The magnitude of the effect depends on (1) the 
amount of H 2 O present in the initial mixture and the corresponding amounts of 
free radicals such as OH present in an equilibrium free stream, and (2) the 
third-body efficiency of H 2 O through its effect on the rate of reaction of 
H + O^ + M ->• HO^ + M. Graphical representations of these two influences are 
presented in figures 8 and 9, respectively. In figure 8, the sensitivity of 
Tig to mass fractions of H 2 O of 0 (air), 0.145, and 0.245 in an initial mix- 
ture at Tq = 1000 K, is presented for conditions of EFS (dashed lines) and 
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with H 2 O only (solid lines). 


At a pressure of 0.2 atm the calculated ignition 
times for all cases range over ±20 percent; however, as the pressure is 
increased to 1.0 atm, the presence of H 2 O alone tends to increase by a 

" ■ " . ^ _ 0.145 and up to 30 for ^ O ” 0.245. 


factor of up to 10 for 


This 


increase in is due to the depletion of H atoms by the previously men- 

tioned hydroperoxyl (HO 2 ) chain terminating reaction. When equilibrium free 
stream is used, the ignition times are reduced as pressure increases to 
1.0 atm by a factor of about 3, relative to the case with H 2 O only. This 
reduction is due to the increased levels of OH. 


The sensitivity of ignition times to the third-body efficiency of H 2 O, 
i.e., m^ ^ - in the reaction H + 0_ + M > HO^ + M at = 1000 K is pre- 

2 2 C 


sented in figure 9. The con^uter program currently uses a value ^04” 


relative to the third-body efficiency of N 2 * A change in ”^04 

preceding reaction causes a corresponding change in the rate of the reaction 
with H 2 O as the third body. The effect of increasing and decreasing m^^ q ^ 
by a factor of 2 is illustrated in figure 9 for vitiated air (EFS) with ^ 

Q = 0.145 at Tq = 1000 K. Also shown is the case for air (EFS) with no 


H 2 O. At pressures near 0.2 atm, the presence of H 2 O with its high third-body 
efficiency has little effect on values. As the initial pressure 

approaches 1.0 atm, the calculated ignition times show a large sensitivity to 
the value of O 4* doubling ^04 increases by a factor of 10 at 

1 atm; halving ^04 ^®^nces by a factor of 2 and gives values of 


T. near those of H 2 "air (EFS). These results emphasize the importance of 
knowing the efficiency of H 2 O relative to N 2 in modeling the chemistry and in 
predicting ignition times of H 2 in a vitiated air stream at static tempera- 
tures of 1000 K or less and pressures near 1.0 atm. 


The sensitivity of ignition times to varying concentrations of NO is pre- 
sented in figure 10 for H 2 "air and H 2 -vitiated air at T^ = 910 K. For H 2 “air 
in figure 10(a) at low pressures (near 0.2 atm), a concentration of NO of as 
much as 0.01 by mass has little effect on the ignition time, causing a reduc- 
tion of less than 20 percent in for pure air. As the initial pressure 

approaches 1.0 atm, however, an NO mass fraction of 0.01 reduces by a 

factor of 50. For an NO mass fraction of 10"^ (10 ppmw), which is the NO 
concentration in air in equilibrium at 910 K, is reduced by about a 

factor of 4. The experimental results reported in reference 15 indicated 
reductions in by factors of 2/3 and 1/2 at pressures of 0.27 and 

0.5 atm, respectively, due to molar concentrations of NO^ between 0.005 and 
0.01. For a stoichiometric H 2 “air mixture, the mass fraction of NO^^ is about 
1.5 times the mole fraction. These results, indicated by squares in fig- 
ure 10(a), agree reasonably well with the present calculations. Data for 
shock-tube ignition of H 2 “air at 910 K (ref. 36), denoted by the circles in 
figure 10(a), also show good agreement with the preset calculations. 
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Figure 10(b) shows that for H 2 ~vitiated air, increasing the mass fraction 
of NO in the initial mixture from 10“^ (the EFS value) to 0.01 causes signifi- 
cant reductions in ^ig' particularly at pressures greater than 0.5 atm. For 
an NO mass fraction of 0.01, the ignition times for both H 2 ~air and H 2 - 
vitiated air are nearly equal. The chemical mechanism responsible in part for 
these reductions is the reaction NO + HO 2 $ OH + NO 2 , which produces the 
active chain carrier OH through the depletion of HO 2 . Because HO 2 is formed 
through a chain terminating reaction that depletes H atoms, it is responsible 
for the retarding of ignition in vitiated-air mixtures. A more complete mech- 
anism for the effect of NO^ on H 2 — air ignition is proposed in reference 37 . 

The results in figure 10 indicate at strong influence of NO on ignition times 
at low temperature and pressures above 0.5 atm. 


Effect of other chemical additives .- Because of the difficulty of obtain- 
ing ignition in a short time (or distance) at the low temperatures (850 to 
1000 K) and pressures (0.30 to 1.0 atm) typical of a scramjet combustor, cer- 
tain chemicals have been considered to prompt ignition. These chemicals, such 
as H 2 O 2 and O 3 , would supply OH free radicals and active oxygen atoms which 
would shorten the ignition time. The oxides of nitrogen or other nitrogen 
compounds which produce these molecules could be used over a limited range of 
pressures and temperatures. The reduction in ignition times for stoichiomet- 
ric mixtures of H 2 — air and H 2 — vitiated air at 1.0 atm obtained by adding mass 
fractions of O^, H 2 O 2 , NO, and NO 2 of 0.01 is presented in figures 11(a) 
to 11 (d), respectively. The reduction in x. is given by the parameter X 

defined as the ratio of x^^g with the chemical additive to Xj^g from fig- 
ure 6 . Chemicals and H 2 O 2 provide substantial reduction in ignition times 
at all temperatures for both air and vitiated air. Although not presented, 
reduction of the pressure to 0.5 atm was observed to have little effect on the 
reduction of x^ for both O 3 and H 2 O 2 . For the NO^^ additives in fig- 
ures 11 (c) and 11 (d), the significant reduction occurs only at the lower 
temperatures (less than 1100 K) . This large reduction in x. at the low 
temperatures occurs because of the same chemical mechanism discussed in the 
previous section - the formation of OH through the depletion of HO, reacting 
with NO. 


The effect of atmospheric O 3 encountered in the flight of a scramjet 
engine also was examined. At altitudes of about 30 km the O 3 concentration is 
approximately 10 ppnw. Peak O 3 concentration is about 13.5 ppmw at 36 km. 

The reduction in ignition time for a trace O 3 concentration of 10 ppmw at 
1 atm is given in figure 12. Substantial reduction occurs at the lower ini- 
tial temperatures (less than 1000 K). Lowering the pressure from 1.0 to 
0.5 atm did not change the result. 


Reaction Time 

In the present analysis the reaction time is defined as the time differ- 
ence between the ignition time and the total reaction time required to achieve 
95 percent of the equilibrium temperature rise (see eqs. (30) and (31)). The 
effect of pressure on the variation of reaction time with temperature is pre- 
sented in figure 13 for stoichiometric mixtures of H 2 -air and H 2 ~vitiated 
air. In contrast to the ignition time, vitiated air has a negligible effect 
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on the reaction time. In the semilog plot of figure 13, the variation of 
with T is linear and shows a strong dependence on p. Ihe present calcula- 
tions also indicate that is independent of ( 1 ) the initial concentration 

of free radicals and atoms (similar results were reported in refs. 7 and 8 )? 

(2) the presence of chemical additives such as NO, ^2^2^ snd and 

( 3 ) within ±10 percent, the mixture equivalence ratio in the range 0.5 

to 1.5. These reaction times can be correlated, within ±10 percent, by the 
equation 


T 

R 


325p 


■ 1.6 


exp(-0.8T^/1000) 


(34) 


A comparison of the correlation given by equation (34) and reaction times 
from references 4, 7, and 8 is given in figure 14 for a stoichiometric H 2 "air 
mixture at p = 1.0 atm. The large difference between the present results and 
the earlier work of reference 4, which gave reaction times as much as 5 times 
faster than the present theory, can be traced back to the rate constants for 
three— body recombination exothermic reactions. 

Effect of changes in reaction-rate constants .- One reaction which has a 
major effect on the calculated reaction time is H + OH + M + H 2 O + M. The 
reaction rate constant k 3 for this reaction in the present chemical model is 
similar to the value recommended in references 20 and 21 with N 2 as the third 
body with an efficiency of 1. However, the rate data of references 20 and 21 
exhibited considerable scatter, with a factor of 7 representative of the upper 
limit change in the rate constant. This factor of 7 was used to increase the 
rate constant (k^) and reaction times recomputed. No estimate of the lower 
limit change was attempted. The comparative effect on of this change xn 

k^ by a factor of 7 is presented in figure 15 for stoichxometric mixtures of 
H -air and H 2 ~vitiated air at pressures of 0.5 and 1.0 atm. For each pres- 
sure, the upper curves are for the value of k^ given in table I; the lower 
curves are for a rate constant of 7 times k^* both cases, the third— body 

efficiency of H^O was n\H^O,3 "" ^ table I). The effect of the factor of 

7 and the high efficiency of H 2 O mean that the rate constant for^the reaction 
H + OH + H 2 O > H 2 O + H 2 O will have a coefficient = 2.2 x 10 compared 

with the table I^^value of A 3 = 0.52 x for the reaction 

H + OH + N 2 ^ H 2 O t N 2 « result, shown in figure 15, is a reduction in 

reaction time by a factor of about 3. The calculated values of with the 

faster rate constants agree better with previous results and fall between the 

T values from references 4 and 7 (compare fig. 14 with fig. 15(b)). 

R 

Effect of third-body efficiency of H 2 O.- Another way to affect reaction 
time is to change the reaction rate through the third-body efficiency of H 2 O 
in the reaction H + OH + H 2 O ^ H 2 O + H 2 O. In reference 38 this reaction is 
considered the dominant path for the production of H 2 O, with a rate constant 
of approximately 2.7 x 10 cm^-mol“^-s ^ at 1650 K. In reference 39, the 
rate constant had to be 2.4 x 10 ^^ exp( 500/RT )cm^-mol’^-s to agree with 

their data. In the present chemical mechanism the rate constant, ^including^a 
third-body efficiency factor of 6 for H 2 O, is2.6xio T cm -mol -s • 
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This value is approximately the value recommended by references 20 and 21. A 
comparison of the three rate constants for the production of H-O is given in 
the following table: ^ 


To, K 

Source 

8 — ^ 1 

cm -mol -s” 

1500 

Ref. 38 

2.7 X 10 ''^ 

1 

Ref. 39 

to 

. 

00 

X 

0 

1 

Present study 

4.5 X 10 '^^ 


In order to check the sensitivity of to a change in this rate of reaction, 

the third-body efficiency of H 2 O was varied from 1 to 40, with the latter value 
corresponding to the rate constants from references 38 and 39. Note that the 
value of the rate constant for the present calculations includes the factor of 
6 for the third-body efficiency of H 2 O. The effect on Tj^ of varying 
"* 820,3 I to 40 is presented in figure 16 for T^ = 1500 K, p = 1 atm, 

and stoichiometric mixtures of H 2 -air and H 2 -vitiated air. These results are 
represented by the parameter F, defined as the ratio of T_ to x for 

fxgure 13, with "^^ 0,3 ^ spread between the H 2 -air and 82 - 

vitiated air reaction times disappears when the third-body efficiencies of 8 O 
and N 2 are equal 3 = 1 j and increases to a difference of a factor of 2 ^ 

when rtijj^o,3 "" giving a rate constant of 3.0 x 10 ’'^ cm^-mol~^-s“^ . 

It is apparent from the comparisons of figures 15 and 16 that the rate con- 
stants for the recombination of 8 and OH to form H 2 O will have an important 
impact on reaction times in scram jet combustors. These results point out the 

need for more accurate rate data for this reaction with N, and 8,0 as third 
bodies. ^ ^ 


CONCLUDING REMARKS 

to analytical study of the chemical mechanism involved in the ignition and 
reaction of the 82 -air system has been performed at conditions representative 
of a scramjet combustor. The calculations were made with a one-dimensional 
kinetics program with an H 2 -O 2 -N 2 finite-rate chemistry model consisting of 
60 reactions in 20 species. An objective of these calculations was to provide 
information for an assessment of the conditions at which kinetics are important 
in scramjet combustors. Specifically, the study considered how ignition and 
reaction times are affected by initial temperature and pressvire; fuel equiva- 
lence ratio; the presence of contaminants such as 82 ©, 8 , 08, and NO ; chemical 
additives H 2 O 2 and O 3 ; and the reaction rate of key reactions. ^ 
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Results of the calculations indicated that, within a 20 percent error, 
ignition time is nearly constant over a range of equivalence ratios from 0.5 
to 1.7. Relative to air in chemical equilibrium, the presence of H 2 O in the 
initial mixture (as from a vitiated main streeun) increased ignition times at 
initial temperatures below 1000 K and reduced ignition times by about 33 per- 
cent at temperatures above 1200 K. These results were found to be due to the 
depletion of H atoms by the accelerated production of HO 2 at low initial tem- 
peratures through the reaction H + O 2 + H 2 O ^ HO 2 + H 2 O, because the third-body 
efficiency of H 2 O may be 13 times that of N 2 . At higher ten^eratures the dis- 
sociation of H 2 O produces OH concentrations that counteract the H atom deple- 
tion and reduce ignition time. The presence of chemicals such as H 2 O 2 , 0^, NO, 
and NO 2 , in concentrations (mass fractions) of 0.01 in both H 2 ~air and H 2 - 
vitiated air mixtures, significantly reduced ignition times at low temperatures 
(near 910 K) . At high temperatures (above 1100 K) NO and NO^ had little effect 
on ignition times, whereas H 2 O 2 and O 3 still caused a reduction of ignition 
times by about a factor of 7. 

The present calculations indicate that reaction time, defined as the time 
from the onset of ignition to attainment of 95 percent of the equilibrium tem- 
perature, is independent of the concentration of free radicals, chemical addi- 
tives, and within 10 percent, the equivalence ratio between 0.5 and 1.5. In 
addition, the presence of H 2 O in a vitiated main stream at mass fractions of up 
to 0.245 had very little effect on reaction time. However, a conparison of 
rate constants for the reaction H + OH + M + H 2 O + M indicates a need for more 
accurate rate data with N 2 and H 2 O as third bodies. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
June 3, 1981 


21 


REFERENCES 


1. Guy, Robert W.; and Mackley, Ernest A.: Initial Wind Tunnel Tests at 

Mach 4 and 7 of a Hydrogen-Burning, Airframe-Integrated Scramjet. NASA 
presented at the 4th International Symposium on Air Breathing 
Engines (Lake Buena Vista, Fla.), Apr. 1-6, 1979. 

2. Bittker, David A.; and Scullin, Vincent J. ; General Chemical Kinetics 

Computer Program for Static and Flow Reactions, With Application to Com- 
bustion and Shock-Tube Kinetics. NASA TN D-6586, 1972. 

3. McLain, Allen G. ; and Rao, C. S. R.: A Hybrid Computer Program for Rapidly 

Solving Flowing or Static Chemical Kinetic Problems Involving Many Chemi- 
cal Species. NASA TM X-3403, 1976. 

Harold S. : A Theoretical Analysis of Non— Equilibrium Hydrogen- 

Air Reactions in Flow Systems. [Preprint] 63113, American Inst. 

Aeronaut. & Astronaut., Apr. 1963. 

5. Jachimowski, Casimir J. ; and Houghton, William M. : Effect of Carbon Diox- 

and Water Vapor on the Induction Period of the Hydrogen— Oxygen Reac- 
tion. NASA TN D-4685, 1968. 

6. Brokaw, Richard S. : Analytic Solutions to the Ignition Kinetics of the 

Hydrogen-Oxygen Reaction. NASA TN D-2542, 1964. 

7. Erickson, Wayne D. ; and Klich, George F.: Analytical Chemical Kinetic 

Study of the Effect of Carbon Dioxide and Water Vapor on Hydrogen- Air 
Constant-Pressure Combustion. NASA TN D-5768, 1970. 

8. Carson, George T., Jr.: Analytical Chemical Kinetic Investigation of the 

Effects of Oxygen, Hydrogen, and Hydroxyl Radicals on Hydrogen-Air Com- 
bustion. NASA TN D-7769, 1974. 

9. Bascombe, K. N. : Calculation of Ignition Delays in the Hydrogen-Air Sys- 

tem. Combust. & Flame, vol. 11, no. 1, Feb. 1967, pp. 2-10. 

10. Skinner, Gordon B.; and Ringrose, Gordon H. : Ignition Delays of a 

Hydrogen-Oxygen-Argon Mixture at Relatively Low Temperatures. J. Chem. 
Phys., vol. 42, no. 6, Mar. 15, 1965, pp. 2190-2192. 

11. Abbett, M. A.: The Effects of Water on the Combustion History of a 

Hydrogen-Air System. Tech. Rep. 427, Gen. Appl. Sci. Lab., Inc., Mar. 
1964. 

12. Rhodes, R. P. ; The Effect of Non-Equilibrium Free-Radical Concentration on 

Ignition Delay in the Hydrogen-Air System. AEDC-TDR-64-24 1 , U.S. Air 
Force, Nov. 1964. (Available from DTIC as AD 608 566.) 

13. Edelman, Raymond B.; and Spadaccini, Louis J. ; Analytical Investigation of 

the Effects of Vitiated Air Contamination on Combustion and Hypersonic 
Airbreathing Engine Ground Tests. AIAA Paper No. 69-338, Apr. 1969. 


22 



14. Bahn, Gilbert S. : Calculations on the Autoignition of Mixtures of Hydrogen 

and Air, NASA CR-1 12067/ 1972, 

15. Slack, M,; and Grille, A,: Investigation of Hydrogen-Air Ignition Sensi- 

tized by Nitric Oxide and by Nitrogen Dioxide, NASA CR-2896, 1977. 

16. Engineering Staff, AiResearch Manufacturing Co.: Hypersonic Research 

Engine Project - Phase II. Aerothermodynamic Integration Model Develop- 
ment. Final Technical Data Report. NASA CR-132654, 1975. 

17. Momtchiloff, I. N.; Taback, E. D.; and Buswell, R. F.: Kinetics in 

Hydrogen— Air Flow Systems. I. Calculation of Ignition Delays for 
Hypersonic Ramjets. Ninth Symposium (International) on Combustion, 
Academic Press, Inc., 1963, pp. 220-230. 

18. Schott, G. L.; and Kinsey, J. L. : Kinetic Studies of Hydroxyl Radicals in 

Shock Waves. II. Induction Times in Hydrogen-Oxygen Reaction. J. Chem. 
Phys., vol. 29, no. 5, Nov. 1958, pp. 1177-1182. 

19. Jensen, R. M.; Bryce, C. A.; and Reese, B. A.: Nonequilibrium Combustion 

Effects in Supersonic Streams. NASA CR-120932, 1972. 

20. Baulch, D. L. ; Drysdale, D. D. ; Horne, D. G. ; and Lloyd, A. C.: Evaluated 

Kinetic Data for High Temperature Reactions. Volume 1 - Homogeneous Gas 
Phase Reactions of the H 2 -O 2 System. CRC Press, c.1972. 

21. Baulch, D. L.; Drysdale, D. D.; and Horne, D. G. : Evaluated Kinetic Data 

for High Temperature Reactions. Volume 2 - Homogeneous Gas Phase Reac- 
tions of the H 2 -N 2 -O 2 System. Butterworth & Co. (Publ.), Ltd., c.1973. 

22. Hindmarsh, A. C. : GEAR: Ordinary Differential Equation System Solver. 

UCID-30001, Rev. 1, Computer Documentation, Lawrence Livermore Lab., 

Univ. California, Aug. 20, 1972. 

23. Gear, C. William: Numerical Initial Value Problems in Ordinary Differental 

Equations. Prentice-Hall, Inc., c.1971. 

24. Wakelyn, N. T. ; and McLain, Allen G. : Polynomial Coefficients of Thermo- 

chemical Data for the C-H-O-N System. NASA TM X-72657, 1975. 

25. Slack, M. W. : Rate Coefficient for H + O 2 + M = HO 2 + M Evaluated From 

Shock Tube Measurements of Induction Times. Combust. Flame, vol. 28, 
no. 3, 1977, pp. 241-249. 

26. Jensen, D. E.; and Jones, G. A.: Reaction Rate Coefficients for Flame 

Calculations. Combust. & Flame, vol. 32, no. 1, May 1978, pp. 1-34. 

27. Jachimowski, Casimir J. ; and Houghton, William M. : Shock-Tube Study of the 

Induction-Period Kinetics of the Hydrogen-Oxygen Reaction. Combust. & 
Flame, vol. 15, no. 2, Oct. 1970, pp. 125-132. 


23 


28. Brokaw, Richard S. : Rate of Reaction Between Molecular Hydrogen and 

Molecular Oxygen. NASA TM X-2707, 1973. 

29. Dixon-Lewis, G. ; Greenberg, J. B.; and Goldsworthy, F. A.: Reactions in 

the Recombination Region of Hydrogen and Lean Hydrocarbon Flames. Fif- 
teenth Symposium (International) on Combustion, Combustion Inst., c. 1975 
pp. 717-730. 

30. Getzinger, R. W.; and Blair, L. S.; Recombination in the Hydrogen-Oxygen 

Reaction; A Shock Tube Study With Nitrogen and Water Vapour as Third 
Bodies. Combustion & Flame, vol. 13, no. 3, June 1969, pp. 271-284. 

31. Halstead, C. J.; and Jenkins, D. R. ; Rates ofH+H+M+H +M and 

H + OH + M ->• H 2 O + M Reactions in Flames. Combust. & Flame, vol. 14, 
no. 3, June 1970, pp. 321-323. 

32. Lloyd, A. C.: Evaluated and Estimated Kinetic Data for Phase Reactions of 

the Hydroperoxyl Radical. Int. J. Chem. Kinet., vol. 6, no. 2, Mar. 

1974, pp. 169-228. 

33. Russin, William Roger; Performance of a Hydrogen Burner To Simulate Air 

Entering Scramjet Combustors. NASA TN D-7567, 1974. 

34. U.S. Standard Atmosphere, 1976. NOAA, NASA, and U.S. Air Force, Oct. 1976. 

35. Evans, J. S.; and Schexnayder, C. J., Jr.: Critical Influence of Finite 

Rate Chemistry and Unmixedness on Ignition and Combustion of Supersonic 
H 2 -Air Streams. AIAA Paper No. 79-0355, Jan. 1979. 

36. Snyder, A. D. ; Robertson, J. ; Zanders, D. L. ; and Skinner, G. B: Shock 

Tube Studies of Fuel-Air Ignition Characteristics. AFAPL-TR-65-93, U.S. 
Air Force, Aug. 1965. (Available from DTIC as AD 470 239.) 

37. Ashmore, P. G.; and Tyler, B. J.: The Nature and Cause of Ignition of 

Hydrogen and Oxygen Sensitizied by Nitrogen Dioxide. Ninth Symposium 
(International) on Combustion, Academic Press, Inc., 1963, pp. 201-209. 

38. Bulewicz, E. M.; and Sugden, T. M. ; The Recombination of Hydrogen Atoms 

and Hydroxyl Radicals in Hydrogen Flame Gases. Trans. Faraday Soc., vol. 
54, pt. 12, no. 432, Dec. 1958, pp. 1855-1860. 

39. Jenkins, D. R.; Yumlu, V. S.; and Spalding, D. B. ; Combustion of Hydrogen 

and Oxygen in a Steady-Flow Adiabatic Stirred Reactor. Eleventh Sympo- 
sium (International) on Combustion, Combustion Inst., 1967, pp. 779-790. 


24 



TABLE I.- HYDROGEN-AIR-REACTION MODEL WITH THIRD-BODY EFFICIENCIES 
[Rate constant is given by k = AT^ exp ( -Activation energy/ 1 . 987T ) ] 


REACTION REACTION 

NUHBFR 


1 

M 

t 

1*02 


1*0 

4 

? 

H 

♦ 

1*HZ 

= 

1*H 

4 

T 

M 

♦ 

1*H2Q 

= 

1*H 

4 

U 

1*H 

+ 

1*02 

= 

1*H02 

4 

5 

M 

«- 

1*N02 

= 

1*N0 

4 

6 

M 

♦ 

l*N0 

- 

l*N 

4- 

7 

1*0 


1*C0 

= 

1*CQ2 

4 


1*H 

+ 

1*N0 

= 

1*HN0 

4 

9 

M 

♦ 

1*H202 

= 

1*0H 

4 

10 

1 *0H 


1 *N0 

* 

1*HN02 

4 

11 

1 *0H 

♦ 

1*N02 

= 

1*HN03 

4 

1? 

H 


1*03 

= 

1*02 

4 

13 

M 

♦ 

l*HC0 

= 

1*00 

4 

14 

1*0 

♦ 

1 *H 

= 

1*0M 

4 

15 
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♦ 

1*0 

= 

1*0H 

4 

16 

1*H2 

+ 

1*0H 

» 

1*H20 

4 
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1*H 
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1*CH 

4 
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♦ 
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(a) H2’-air; p = 0,5 atm 


Figure 1.- Species mass fraction as function of time, 4^ = 1*0; T = 1000 K 













Figure 4. - Ignition times for 















(b) Effect of temperature 


Figure 6.- Concluded. 






P ♦ atm 

Figure 7.- Effect of equilibrium free strecim on ignition times 

T = 1300 K; 4> = 1.0. 
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Figure 16.- Effect of third-body efficiency of H 2 O in reaction 
H+OH+M-^HO + Mon reaction time. = 1500 K; 

n = 1 atm; 4> = 1.0. 
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